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ABSTRACT 


In this paper, a robust mathematical model was developed to optimize 
production, distribution, and inventory planning for a perishable supply chain 
consisting of three factories, one warehouse, and three retailers to maximize 
profit. The proposed Fig. Imodel was formulated using mixed-integer linear 
programming and solved using Xpress[VE. The model considers the 
disruption of factory capacity due to uncertain interruptions. In addition, it 
deals with a sustainable supply chain network and interruptions that affect 
factory capacity while maintaining the FIFO inventory policy to maximize the 
network’s overall profit and reduce the amount of expired products to reduce 
the environmental impact. The model successfully handles the problem of 
designing and planning a long-term supply chain network for perishable 
products. The capability and accuracy of the model were verified by solving 
several scenarios. The effect of the minimum and maximum expected values 
of the interruption was studied. It may be concluded that the robust SCN 
design is affected mainly by the value of the maximum expected interruption. 
The design and performance of the SCN are affected by the capacity 
disruption of the factories, and the model succeeded in designing and 
planning the SCN under all conditions. 


Keywords: Supply Chain; Perishable; Robust; MILP; XpressIVE 


1. INTRODUCTION 
Sustainable supply chain design (SSCD) has three pillars: economic, 


environmental, and social. Organizations are now attempting to strike a 
balance between the triple bottom line of sustainability (TBLs) [1], [2]. The 
term “sustainable” was coined at the end of the 1980s, and the term 
“development” was coined to address economic, social, and environmental 
concerns [3]. In today’s world, sustainability is a major factor in a variety of 
businesses. However, few mathematical models combine the economic and 
social consequences of environmental factors and the supply chain [1]. 

For perishable goods, Jiang et al. [4] considered the carbon footprint when 
designing the SSCN. They created a bi-objective mathematical model to 
determine the best partner selection strategy, technology, and transportation 
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method, as well as raw material and recovery material supply. In the sphere of the global food system, Rohmer et al. [5] offered a 
design approach that included both cost and environmental objectives. Yakavenka et al. [6] created a multi-objective model for 
perishable food SSCD. Their suggested model helps in the selection of transportation modes, routes, facility sites, and related 
product flows. Eskandari Khanghahi et al. [7] created a model of multi-objective functions for the blood supply chain that took 
sustainability into account. 

To choose acceptable suppliers in the field of pharmaceutics, Halim et al. [8] and Tavakkoli Moghaddam et al. [9] established a 
mathematical model for the reverse supply chain of perishable foodstuffs, which they used for a meat product manufacturing 
company. Varsei and Politkovskay [10] created a multi-objective model to design a long-term wine supply chain. Zahiri et al. [11] 
proposed a multi-objective sustainable model for constructing a pharmaceutical supply chain structure that takes strategic and 
tactical actions into account. 

Alashhab and Mlybari [12] developed a robust green supply chain network planning optimization model that considers possible 
risks to identify the production, inventory, and shipping methods. Jouzdani et al. [13] developed a dynamic SCND that takes into 
account traffic congestion and demand uncertainty. They used a case study to demonstrate the applicability of their concept in the 
dairy business. Pishvaee et al. [14] presented a case study of an accelerated Benders decomposition technique for sustainable supply 
chain network design in an uncertain supply chain for medical needles and syringes. 

For the perishable food supply from a single manufacturer with a single warehouse, a proposed solution was to apply Monte 
Carlo simulation to the perishable inventory routing problem (PIRP) [15]. The sustainable supply chain problem regarding the 
behavior of demand and the number of locations was examined by Al-Ashhab et al. [16], who determined that the final objective 
was to keep up with the demands of the retailer while using minimum travel between the warehouse and the retailer. To solve the 
supply chain problem, they used the nonlinear generalized reduced gradient solver, which considers the perishable product in 
multi-periods and prevents the expiry of the product. 

Several studies have shown that hierarchical production planning is a useful strategy to deal with the difficulties associated with 
multilevel production planning and scheduling issues that arise in real-world systems [1]. Yakavenka et al. [6], and [17] have used 
multi-criteria decision-making (MCDM) techniques. Repoussis and Gounaris [18], and Ferguson and Koenigsberg [19] have 
addressed a variety of assumptions about the industrial environment. When it comes to production planning in the face of 
uncertainty, Thiel [20] offers an in-depth look at mathematical models to do so. 

Mohammadi et al. [21] addressed the dynamic multilevel lot size issue in a broad product design with single and multiple 
limiting resources, as well as time frames. Kim and Bell [22] looked at the difficulties associated with planning multi-item, 
multilevel backlog production and found them to be very difficult. Lim and Shanthi Kumar [23] developed a multilevel lot size 
difficulty model with backlogs to simulate capacity-constrained multilevel lot size problems with backlogs. Chan [24] investigated 


multi-product, multi-period, and multi-stage production planning using sequence-dependent settings in a multi-product, multi- 


period environment. 


Fig. 1 The general structure of the proposed supply 


In this research, the authors have decided to develop a robust supply chain network model considering perishable products using 
mixed-integer linear programming that integrates both financial and environmental measures for designing and planning a supply 
chain network. This is carried out in consideration of the uncertain interruption of machines. The main objective of the model is to 
maximize profit by maximizing the total revenue and minimizing the total cost (including expiration costs) to prevent the expiration 


of perishable products and to protect the environment. The proposed supply chain network is shown in Fig. 1. 
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The following assumptions have been considered: 

The supply chain is a single product. 

The retailer’s demand is known. 

The shelf life of products is known and limited. 

The first products produced are issued first, according to FIFO inventory management. 

There is no initial inventory. 

The expired products are considered waste and disposed of. 

The inventory is kept only in the warehouse. 

Backlog is allowed, and the unsatisfied demand through one period can be satisfied in any of the following periods. 


2. RESEARCH SIGNIFICANCE 


The proposed model helps supply chain managers to design and plan a robust SCN of the perishable products to assure 
sustainability by maximizing the total profit and avoiding the expiration of products to save environment taking into account the 
uncertainty of production capacity due to the uncertain interruption in the working environment. The work done in this scope is 


very rare as presented in the literature review. 


3. MODEL FORMULATION 


The indices, parameters, and variables used in the proposed model are: 


Indices 
t Index of planning periods (t = 1, ..., T) 
a Index of the ages of products (a= 1, ..., A) 
f Index of potential factories (f = 1, ..., F) 
r Index of potential retailers (r= 1, ..., R) 
Parameters 
SL The product's shelf-life 
Drt The demand of a retailer (r) in period (t) 
CapMmt The capacity of a factory (f) in period (t) 
CapW, The capacity of a warehouse in period (t) 
Upper The maximum expected value of interruption 
Lower The minimum expected value of interruption 
pr Price per unit 
TC Transportation cost per unit per distance 
DMW,, Distance between a factory (f) and a warehouse 
DWC, Distance between a warehouse and a retailer (c) 
IHC Inventory holding cost per unit per period 
EC Expired cost per unit 
SC Shortage cost per unit 
MC Material cost per unit 
PC Production cost per unit 
a Fixed cost 
Variables 
QFW, Quantity of units transferred from a factory (f) to a warehouse in a period (t) 
QW Cet Quantity of units transferred from a warehouse to a retailer (r) in a period (t) 
lat Inventory of the age (a) at the end of a period (t) 
Texpt Expired inventory at the end of a period (t) 
Interruption The percentage of the outage time to the total planned time (uncertain) 
Xf The binary variable equals 1 if a factory (f) is established; 0 otherwise 
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Z Total profit 


3.1. Objective Function 


The profit of the supply chain network is maximized using the economic goal function provided in Eq (1). 


Profit = Total revenue—Total costs 
Total costs include fixed, production, material, transportation, holding, expiration, and shortage costs 


The suggested model’s objective function is as follows: 


Z= y > QWR,, pr — >, FCX, — yy >» QFW, (PC + MC) — >. >. QFW,, TC DFW; — >» > QWR,, TC DRW, 
r ot £ m t f t r ot 
t t 
- YY wath = Yi texns tc - YY se(’n.- Yawn.) 
a ot t r t 1 1 


(1) 


The objective function may be broken down as follows: 
The sales revenue is calculated using Eq. (2) by multiplying the number of units transferred from the warehouse to the retailers 
in all periods by the unit price. 
(2) 
Sales revenue = >. yy QWR,-, pr 
© t 
The fixed cost of the factory is calculated using Eq. (3) by multiplying the fixed cost per factory by a binary decision variable Xf 


that equals 1 if a factory is established at the candidate location f and 0 otherwise, as presented in Eq. (4). 


Fixed cost of the factories = y FC X¢ (3) 
f 


1 If a factory (f) is established 
Xp -| (4) 


0 otherwise 


The production cost is equal to the number of units transferred from the factory to the retailers (r) during all periods (t) 


multiplied by the production cost per unit, as presented in Eq. (5) 
Production cost = » > QFW; PC (5) 
f ¢ 


The cost of the material is equal to the number of units transferred from the factory to the warehouse during all periods 


multiplied by the cost of the material per unit, as presented in Eq. (6). 
Material cost = oe y QFW;, MC (6) 
ft 


The transportation cost is calculated using Eq. (7) by multiplying the number of units moved from each facility by the 


transportation cost per unit distance and the distance between them. 
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Transportation cost = > > QFW,;, * TC DFW; + > y QWR,, TC DRW, (7) 
iA t r t 


Multiply the remaining inventory in all periods by the unit inventory holding cost to obtain the inventory holding cost, as 


presented in Eq. (8). 
Inventory holding cost = >, y Tq [HC (8) 
a t 


The expiration cost is calculated using Eq. (9) by multiplying the expired inventory at the end of all periods by the unit 


expiration cost. 
Expiration cost = oe Texpsy EC (9) 
t 


The difference between the cumulative demand and the cumulative received amounts of each retailer is multiplied by the unit 


shortage cost per period to compute the shortage cost as presented in Eq. (10). 
t t 

Shortage cost = ». ». SE (> Di - y awn) (10) 
r ot 1 1 


[1] Constraints 
The model constraints are presented in this section. 


Balancing Constraints 


lot = y QFWp, — y QWR, t=1, a=1 (11) 


fCF reR 
t S=1 
lat = QFW, — max (> awe.) = a Vea | t=2..4ha=1 (12) 
bi a=1 


s-1 
Tat = MAX 4 Ig-4¢-1 — Max >: QWR,;+ | — > Tie-107,07,¢ = 2,..,T,a = 2,...,8 (13) 
j=a 


TreR 
The equivalence of the inventory balance in the warehouse is ensured. The number of units transferred in the first period from 
the factory to the warehouse matches the number of units transferred from the warehouse to the retailer, in addition to the 
inventory as presented in Eq. (11). 
Constraints (12) and (13) are used to determine how much inventory is left in the warehouse after each period. They ensure that 
units are handled by the FIFO inventory policy. Under a FIFO inventory policy, the oldest inventory, or inventory of the highest 
age, is used first to meet demand, followed by the inventory of intermediate age. Ultimately, demand is met by the inventory of the 


freshest units, i.e., inventory of the lowest age. 
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Texpt =OVt< SL (14) 
Texpt = latVt => SL,a= SL (15) 
Constraint (14) determines the inventory that has expired, and it ensures that no product will expire before its shelf life expires. 


Constraint (15) determines that the expired inventory could only appear once the shelf life has been reached. 


Capacity Constraints 
QFW;, S CapFry * (1 — Interruption, ) Xe V ft (16) 


lower < Interruptions, < Upper V f,t (17) 
Constraints (16) and (17) ensure that the number of units transferred from factories does not exceed the capacity of the factories 


calculated and that the interruption does not exceed the upper and lower rates of interruption. 


Y OWRre <Drt,,t=1 (18) 
TER 
t t-1 
OWks= > pee > QWR yy ,t > 2 (19) 
t=1 t=1 


Constraint (18) ensures that the number of units transferred from the warehouse to the retailers during the first period does not 
exceed the amount of demand of the retailers in the first period. 

Constraint (19) ensures that the number of units transferred from the warehouse to the retailers during the periods does not 
exceed the difference between the cumulative demand of all retailers and the cumulative quantity transferred for each period from 


the warehouse to the retailers. 


>. OF Wr <Capw,Vf,t=1 (20) 
Fer 
Ler QFWre + veer lat-1 < CapW,V f,t>1 (21) 


Constraint (20) ensures that the number of units transported from the factory to the warehouse during the first period does not 
exceed the warehouse capacity. 
Constraint (21) ensures that the overall number of items entering the warehouse in each period, as well as products left over 
from the previous period, does not exceed the prespecified capacity of the warehouse. 
QFWt, QWRrts lat, lexp, 20,Vf,1r,a,t (22) 
X, €{0,13,Vf (23) 


Constraint (22) prevents the quantities from taking negative values, and constraint (23) ensures that X;is binary. 
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4. RESULTS AND DISCUSSIONS 


In this section, the effect of the minimum and maximum expected values of interruption on the network performance will be 
studied. In addition, we will study the effect of the capacity disruption on the supply chain network's design, planning, and 
expected profit. The model parameters are assumed as presented in Table . 


Table 1 Working parameters 


Parameter Value 
Shelf life 3 (period) 
CapFrt 6000 (unit) 
DFW, 23 & DWR,123 200, 300, 400 (km) 
CapW, 100,000 Unit 
THC 3 (USD/unit/t) 
EC,MC and PC 10 (USD/unit) 
SC 5 (USD/unit) 
FC 4000 (USD/unit) 
TC USD/unit) 
PR (USD/UNIT) 


4.1. Studying the Effect of the Minimum and Maximum Expected Value of Interruption 

The objective of this experiment is to study the variability range in the effect of the minimum and maximum expected values of 
possible interruptions on the different cost elements. For these test problems, all the parameters are kept fixed, with interruption as 
the only variable. The range of the minimum and maximum expected values of possible interruptions is assumed to be from zero to 
0.8 with a step of 0.2. The resulting incomes and costs of these experiments are presented in Fig. 2. It may be concluded that the 
performance of the network is affected mainly by the maximum value of the expected interruption. 


2,00,00,000 
15276000 

1,50,00,000 
A 10735360 = Min. = 
YN 
5 1,00,00,000 = Min. = 0.2 
sg 6067800 . 
2 = Min. = 0.4 
=  50,00,000 , 
2 1312800 aie = Oe 

P x m Min. = 0.8 


-50,00,000 ; : E -3493360 
Maximum interruption 


Fig. 2. The effect of the minimum and maximum expected interruption on the profit. 
4.2. Studying the Effect of The Capacity Disruption on the SCN Design and planning 


In this section, the effect of the capacity disruption on the design and planning of the SCN is studied. All parameters were kept 
constant as shown in Table 2 except for the maximum percentage of interruption. 
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Table 2 The optimal robust SCN design and plan at different interruption rates for the first scenario. 


2 Produced 
L 7 z v, 
6. & 2 : 
oO =) 
E = 2 5 
5 Y U SG 
= a Fl F2 F3 ia wn 
0 6 6 0 6 0 
0.2 6 4.8 1.2 6 0 
0.4 6 3.6 2.4 6 0 
0.6 6 2.4 2.4 1.2 6 0 
0.8 6 1.2 1.2 1.2 3.6 2.4 


As shown in Table it may be noticed that at zero interruption, the capacity of only one factory is enough to fulfil all retailers’ 
demands. While at the interruption rates of 0.2 and 0.4, the capacity of the factories decreases to 4800 and 3600 units per period 
respectively, so the first factory produces at its full capacity and the second factory opens to satisfy the extra amount of demand. 
There is no need for a third factory. 

At an interruption rate of 0.6 and 0.8, the capacity of the factories decreases to 2400 and 1200 units per period respectively, so the 
first and second factories produce at their full capacity and the third factory opens to satisfy the extra amount of demand. 

From this study, it may be concluded that the interruption ratio has an obvious effect on the design and performance of the 
SCN. 


4.3. Studying the Effect of the Capacity Disruption on the SCN Performance 
In this section, the effect of the capacity disruption on the SCN performance; quantities delivered, stored, and shortages is studied 


at different value of interruption rates. The demand of all retailers at all periods are as shown in Table . 


Table 3 Retailer's demand (in Thousands) 


Period 1 2 3 4 5 6 
Demand 2 6 8 10 12 
Period 7 8 9 10 11 12 
Demand 12 10 8 6 4 2 


As shown in Table , at an interruption rate of 0.2, the capacity of the factories decreases to 14,400 units per period. The demand 
exceeds the total capacity in all periods except the first two periods in which the SCN retains some inventory to satisfy some of the 
extra demand in the next periods, according to the excessive increase in demand compared to the SCN's capacity, the shortage 
quantities increase gradually during the last periods. 


Table 4 The optimal robust plan of the SCN at an interruption rate of 0.2. 


Inventor yy 
“ z z g 7 e & B 
g q 5 a) g s =) 
o 5 pe = 5 o 2 
a fa Rs a al a2 a3 E a rS 
jam 
1 6 14.4 6 8.4 0 0 8.4 
9 12 14.4 12 10.8 0 0 10.8 
3 18 14.4 18 7.2 0 0 7.2 2 
4 24 14.4 21.6 0 0 0 0 2.4 8 
5 30 14.4 14.4 0 0 0 0 18 S 
6 36 14.4 14.4 0 0 0 0 39.6 
7 36 14.4 14.4 0 0 0 0 61.2 
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30 14.4 14.4 0 0 0 0 76.8 

24 14.4 14.4 0 0 0 0 86.4 
10 18 14.4 14.4 0 0 0 0 90 
11 12 14.4 14.4 0 0 0 0 87.6 
12 6 14.4 14.4 0 0 0 0 79.2 


At an interruption rate of 0.6, the capacity of the factories decreases to 7200 units per period. The demand exceeds the total 
capacity in all periods except the first period in which the SCN retains some inventory to satisfy some of the extra demand in the 
next periods, as shown in Table 5. According to the excessive increase in demand compared to the SCN's capacity, the shortage 
quantities increase gradually during the last periods. 


Table 5 The optimal robust plan of the SCN at an interruption rate of 0.6. 


Inventory 


2 5 . a sg Ss # 
a cs a al a2 a3 & = & 
1 6 72 6 12 0 0 12 0 
2 12 72 8.4 0 0 0 0 3.6 
3 18 72 72 0 0 0 0 14.4 
4 24 72 7.2 0 0 0 0 31.2 
5 30 in 72 0 0 0 0 54 
6 36 72 7.2 0 0 0 0 82.8 at 
7 36 72 7.2 0 0 0 0 111.6 a 
8 30 72 7.2 0 0 0 0 134.4 
9 24 72 7.2 0 0 0 0 151.2 
10 18 72 ver, 0 0 0 0 162 
11 12 72 Ta 0 0 0 0 166.8 
2 6 72 yee: 0 0 0 0 165.6 


The effect of interruption on the quantities delivered to retailers is presented in Fig. 3. There is a decrease in the delivered 


quantities due to increases in the number of interruptions. 


=@=— Max interruption 0.0 


35,000 —@=— Max interruption 0.2 
30,000 —=@=— Max interruption 0.4 
2 —=@®=— Max interruption 0.6 
ai == Max interruption 0.8 
20,000 
15,000 
10,000 
9,000 9-9—-0—0—0—0—-0-0-0-0-0-0 
0 Period 


Fig. 3. Effect of interruptions on the quantities delivered to retailers. 
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The effect of interruptions on the residual inventory is presented in Fig. 4. The residual inventories decrease because of the 
increase in interruptions. The effect of interruptions on the shortage quantities is presented in Fig. 5. The shortage quantities 
increase because of the increase in interruptions. 


20,000 > 
z == Max interruption 0.0 
15,000 
a =@=— Max interruption 0.2 
10,000 == Max interruption 0.4 
5,000 =@=— Max interruption 0.6 
0 


Period 


Fig. 4. Effect of interruptions on the residual inventory. 


2,50,000 
=@— Max. = 


200,000 =—@— Max. = 0,2 

=—@— Max. = 0.4 

1,50,000 =@= Max. = 0.6 

2 =O Max. 

1,00,0068 

Oo 
50,000 


0 


Period 


Fig.5. Effect of interruptions on the shortage quantities. 


The Effects of Interruptions on the Supply Chain Network's Profit 


The effect of interruptions on the profit is presented in Fig. 6. The increase in interruptions leads to a decrease in profits in almost all 
cases. 
2,00,00,000 
1,50,00,000 
1,00,00,000 
_30,00,000 


0 


Profi 


0.2 O04 0.6 


-50,00,000 
Max interruptions 


Fig. 6. Effect of interruptions on the profit. 


5. CONCLUSION 

The developed study models a supply chain network consisting of three echelons: number of factories, distribution warehouse, and 
number of retailers for perishable products in a limited multi-period time horizon. The model takes into consideration the 
disruption of factory capacity due to an uncertain interruption. In addition, it deals with a sustainable supply chain network and 


interruptions that affect factory capacity while maintaining the FIFO inventory policy to maximize the network’s overall profit, 
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reduce the amount of expired products, and reduce the impact on the environment. The model successfully handled the problem of 


designing and planning a long-term supply chain network for perishable products while considering inventory in the warehouse 


and the backlog. Moreover, a demand that is not supplied in one period can be met in succeeding periods. This aids in the 


prevention of product expiration by ensuring that only the appropriate quantities are generated with no waste. 


The capability and accuracy of the model were verified by solving several scenarios. The effect of the minimum and maximum 


expected values of interruptions was studied. It may be concluded that robust SCN design is mainly affected by the maximum 


expected interruption. The design and performance of the SCN are affected by the capacity disruption of the factories, and the 


model succeeded in designing and planning the SCN under all conditions. 


For future studies, we recommend considering: 
e multi-products. 
e multi-items. 
e multi-objectives. 
e the reverse logistics of the expired products; and 
e stochastic demand 
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